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O 

o 



C/3 



S. Tsui, Y. Q. Wang, and Y. Y. Xue 
Department of Physics and Texas Center for Superconductivity, 
University of Houston, 202 Houston Science Center, Houston, Texas 77204-5002 

C. W. Chu 

Department of Physics and Texas Center for Superconductivity, 
University of Houston, 202 Houston Science Center, Houston, 
Texas 77204-5002; Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, 
California 94720; and Hong Kong University of Science and Technology, Hong Kong 

(Dated: February 6, 2008) 

The polarity-dependent resistive-switching across metal-Pro. yCao.aMnOs interfaces is investi- 
gated. The data suggest that shallow defects in the interface dominate the switching. Their 
density and fluctuation, therefore, will ultimately limit the device size. While the defects gener- 
ated/annihilated by the pulses and the associated carrier depletion seem to play the major role at 
lower defect density, the defect correlations and their associated hopping ranges appear to dominate 
at higher defect density. Therefore, the switching characteristics, especially the size-scalability, may 
be altered through interface treatments. 
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The renewed interest in various resistive switching 
phenomena is largely driven by recent market demands 
for nano-sized nonvolatile memory devices.^ While the 
current boom of consumer electronics may largely be 
attributed to the successful miniaturization of both 
FLASH chips and mini hard drives, cheaper and smaller 
devices are called for. Various resistive hysteretic 
phenomena are consequently studied with the hope 
that the size limitations associated with the related 
physics/chemistry/technology might be less severe;- Our 
limited knowledge about the mechanisms so far, how- 
ever, makes the evaluation difficult. This is especially 
true for the switching across metal-Pro.yCao.sMnOa 
(PCMO) interfacesi^i'^i^i^ Several models, i.e. bulk 
phase-separation,>3' carrier-trapping in pre-existing metal- 
lic domains,^'* and field-induced lattice defects,*^ have 
been proposed. Each possesses its own distinguishable 
size- limitation, e.g. the statistics of the associated lo- 
cal mesostructures. Here, we report our mechanism in- 
vestigation through both the trapped-carrier distribution 
and their hopping range. Our data suggest that the 
characteristics may largely be engineered through the 
mesostructure of the interfacial defects. 

Bulk PCMO, in great contrast with well known semi- 
conductors, has a rather high nominal carrier concentra- 
tion with its high resistivity mainly attributed to hopping 
barriers^ Local defects, therefore, appear as a natural 
cause of the resistive switching. Following this line of 
reasoning, a domain model has recently attracted much 
attentioniL^ There, a tunneling from the electrode to 
some pre-existing interfacial metallic domains has been 
assumed to be the dominant process. Consequently, the 
carrier-occupation in the domains may change with the 
carrier-trapping during the write pulses, and cause the R- 
switch between an on (low resistance) and an off (high re- 
sistance) state. This is realized through either the change 
of the tunneling probability'' or a doping-induced metal- 
insulator transitioui^ Useful devices based on this mech- 



anism, therefore, should typically be much larger than 
these interfacial domains. It is interesting to note that 
even if the "domains" can be reduced to individual lat- 
tice defects (or small clusters) as in the proposed defect 
modification model,* the fluctuation (inhomogeneity) of 
the defect density still sets a limit for the size scalabil- 
ity just like the dopant fluctuation in Si nano-devices^ 
The defect (domain) density, therefore, requires explo- 
ration, and the interfacial capacitance, C(w), can serve 
to distinguish between these models. 

While the domain model may simulate very diver- 
gent dc I — V characteristics by adjusting the fitting 
parameterSji the measured capacitance is expected to 
be ^-independent with Con = Coff {Con > Coff) for 
tunneling-probability (metal-insulator transition) scenar- 
ios, where Con and Co// are C(w) in the on and off 
states, respectively. In the defect modification model, 
however, the C{uj) measures the net trapped carriers 
in the interface, i.e. with Roff/Ron ~ Coff /Con in 
the space-charge-limited-current (SCLC) region, where 
Ron and Roff are R{lu) in the on and off states, respec- 
tively. It should be noted that the carriers trapped be- 
hind a hopping barrier Vhop respond to a step-disturbance 
as exp[— Wot • exp(— V/iop/fes^)]. Therefore, C(w) cx 

oo 

J kBT[dN/dVhop]/(^duj, i.e. the density of states (DOS) 

at the Fermi level for the defects with the hopping barrier 
lower than Vhop — kBTlii{vQ/u)), where fc^ and vq « 10^^ 
sec~^ are the Boltzmann constant and the trial frequency, 
respectively. An experimental challenge, however, exists 
in separating the interfacial C(w) from the bulk contri- 
bution. We have previously reported a preliminary re- 
sult for an Ag-PCMO interface"* through the traditional 
Cole-Cole procedure, which assumes that all C's and i?'s 
arc cij-independent, under a standard two-leads measure- 
ment configuration. In the SCLC region, the observed 
Coff = 2.2 nF and Con = 1-6 nF clearly contradict 
the domain model but qualitatively agree with the defect 
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FIG. 1: R{i^) (triangle symbols) and C{uj) (square symbols) 
of a Au-PCMO setup at on- (solid) and off- (open) states. 
Bottom left inset: the switching series. Bottom right inset: 
three-lead setup. 



modification model. The data, however, also raise a seri- 
ous concern about the size-limitation. In a sense, the ob- 
served C/qK, 10^^ electrons/cm^, where q — 1.6- 10~^^C 
is the electron charge, is a measure of the trapped carri- 
ers (or the shallow defects near the Fermi level), and puts 
a statistical size-limit on the order of 10-100 nm if the 
switching is due to a change in defect density. There- 
fore, a possible route to further minimizing the device 
size would be to increase the C(a;). A dilemma, how- 
ever, arises: denser shallow defects would also enhance 
the thermal excitation. Across a C{ljj) threshold, the 
SCLC may not be reachable, and the denser defects at 
the off states, functioning as donors, might even enhance 
the conductivity. It is therefore even unclear whether 
samples with a much larger C{uj) are switchable. 

We finally found several metal electrode-PCMO film 
configurations, with the low-w interfacial C > 100 
nF/mm^, a value one order of magnitude higher, al- 
though the conditions for reproducibly synthesizing such 
samples are not yet clear. Repeatable switching has been 
obtained (bottom left inset, Fig.QJ. The Cole-Cole plot, 
however, shows that the Coff/Con is even smaller than 
1, while Roff/Ron > 2, a scenario closer to the domain 
model. To explore the issue, a new procedure is devel- 
oped to directly measure the complex interfacial admit- 
tance 1/R{u!) — iujC{uj). This is done by extending the 
previous three-leads i?-measurement to the ofF-phase part 
through a Solartron SI 1260 impedance/gain-phase ana- 
lyzer (bottom right inset. Fig. Resistors/capacitors 
networks were used to verify that the phase uncertainty 
is less than 0.1°. Both R{lo) and C(a;) of the interface, 
therefore, can be accurately deduced over 10^-10^ Hz. 

The C{uj) and the i?(w) observed undergo a step-like 
jump over a narrow range of 0.1 and 3 MHz, with the 
C-jump occurring at higher lo (Fig. This is simi- 
lar to the Maxwell- Wagner relaxation, -'^^ but differs from 
dielectric Debye relaxation. -'^^ Carrier polarization, i.e. 
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FIG. 2: R{T) between 150 and 300 K. Top (bottom) grey 
symbols: data at off- (on-) states. Lines: VRH fits. Inset: 
the calculated defect distribution in the 0.1 mm^ interface 
layer. Open triangles: the on state. Solid circles: the off 
state. 



trapping and hopping in disordered solidsji^ appears to 
be a natural interpretation. In such models, only the 
hopping paths with all barriers Vhop ^ /cbT In (lj/uo) 
contribute to the apparent conductivity l/i?(w), and the 
Lii)dC{ijj) / duj measures the defect distribution against the 
hopping barriers, dN/dVhop, at the Fermi level. The 
jumps (Fig.0, therefore, suggest a defect mesostructure 
with many short conducting clusters (domains) with the 
intra-cluster Vhop < 0.3 eV separated by slightly higher 
inter-cluster barriers, i.e. around 0.4 eV (inset, Fig.|2J). It 
is interesting to note that the Ron and Roff are distinct, 
i.e. switchable, below 0.1-1 MHz, but the correspond- 
ing dN/ dVhop is practically the same as indicated by the 
Cole-Cole plot. This is very different from the samples 
with smaller C{uj), although both show higher R after 
positive pulses. Also, the same R{uj) above 1 MHz at 
both on and off states indicates a limitation on the read 
speed in future potential applications. 

It should be pointed out that such strong dispersions of 
R{uj) and C{uj) suggest that the domain model is also in 
disagreement with the data. The interfacial hopping in 
the model has been attributed to a single barrier layer be- 
tween the electrode and the proposed interfacial domain, 
and should therefore be cj-independent for lu ^ vq. Thus, 
a more complicated defect network should be invoked to 
accommodate the dispersions. It is also interesting to 
note that the thermo-produced carriers can be directly 
deduced from the C{uj) observed. The calculated value, 
> 10^^ sec~^ below 100 kHz within the interfacial layer, 
is far higher than the reported injected current, < 5 ■ 10^'* 
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sec""'^, at the SCLC region4 This supports the above dis- 
cussion for samples with large C(ijj), and suggests that 
the switching mechanism can be modified through inter- 
face engineering. 

To answer the key question of why the switching can 
still occur, R{T) was measured at both on- and off- 
states (Fig. While the i?(T)'s appear to be almost 
parallel, which makes the switching difficult to under- 
stand if Coff ~ Con, a closer examination shows that 
R = a ■ T ■ exp[— (Tq/T)'''] in the variable-range-hopping 
(VRH) formulation might be a better description, as sug- 
gested by the slight curvatures in Fig.|21 For this particu- 
lar sample, the fitting parameters are 7 — 0.85 (0.6f ) and 
To = 2089 K (8353 K) for the on (off) states, respectively. 
The 7 and Tq are traditionally associated with 1/(1 + d) 
and f /DOS, where d is the dimensionality.^'* The changes 
in both 7 and Tq, therefore, tentatively suggest that the 
arrangement of the defect-structures are changed. The 
switching, in such a case, might be more about the local 
structure than the average defect density. 

To verify this assumption, the small-signal ac R{uj) at 
various dc biases and temperatures was measured (inset, 
Fig. 13) . The normalized I — V characteristics between 
8 and 55 °C are scaled into a single trace as suggested 
by Mott's formula of R/R{E = 0) oc exp{-l ■ E/ksT)^ 
or exp(— Z • E /kBT)^ where I and E are the hopping 
range and the electric field, respectively. The physical 
picture, i.e. the longer the hopping range, the larger the 
field effects will be, is straightforward and independent 
of the hopping details, although the absolute value may 
be affected by the f nm assumed thickness^ of the inter- 
facial layer. Our data, therefore, demonstrate that the 
hopping range is longer at the off state for w < 1 MHz. ft 
is also interesting to note the unusually small / for lo > 1 
MHz (Fig. no . This may appear only if the associated 
hopping barriers are negligible, which is further sugges- 
tive of local defect mesostructure. With high enough 
defect density, the formation of local defect-correlations 
becomes a means to modify the switching performance, 
ft should be noted that, in order to accommodate the ob- 
served Cof f I Con , the conductance may be dominated by 
only a few percolation paths, such that simply enhanc- 
ing C(uj) may not necessarily reduce the size limitations. 
Direct observations of such structures are called for and 
planned. 

In conclusion, we have shown that the interfacial re- 
sistive switch found in metal-PCMO systems is a defect- 
mediated process. Whereas a defect density alteration 
occurs in samples with C{uj < 10^ Hz) < 100 nF/mm^, 
we propose that the change in defect mesostructure may 



cause switching in the samples with larger C{uj). Further 
investigation of the parameters controlling the defect dis- 
tribution will shed light on how best to proceed with the 
nanoscaling and benchmarking of future device applica- 
tions. 
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FIG. 3; Hopping range with respect to frequency for both off 
(open triangles) and on (solid circles) states. The switch is 
determined by low frequency inter-cluster hopping, whereas 
at high frequencies, percolative intra-cluster hopping occurs. 
Inset: raw I — V data at several temperatures show the scaling 
of Mott's formula, where To = 300 K. 
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